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Abstract. The quasineutral limit of compressible Navier-Stokes-Poisson sys- 
tem with heat conductivity and general (ill-prepared) initial data is rigorously 
proved in this paper. It is proved that, as the Debye length tends to zero, the 
solution of the compressible Navier-Stokes-Poisson system converges strongly 
to the strong solution of the incompressible Navier-Stokes equations plus a 
term of fast singular oscillating gradient vector fields. Moreover, if the Debye 
length, the viscosity coefficients and the heat conductivity coefficient indepen- 
dently go to zero, we obtain the incompressible Euler equations. In both cases 
the convergence rates are obtained. 



1. Introduction 

In the present paper we study the quasineutral hmit of compressible Navier- 
Stokes-Poisson system 

9tP + div (pu) = 0, (1.1) 

p{dtu+{u ■ V)u} + VP(p,6') + pV$ = /^Au+ + j/)Vdivu, (1.2) 

cvp{dt9 + (u • V)e} + P{p, 0)diy u = kA9 + z/(divu)2 + 2pM{u) : D(u), (1.3) 

-A^A-J^p-l, (1.4) 

for X e C R^{N — 2, 3), the iV-dimensional torus, where p, u = (iti, . . . , un), 
0, and <I> denote the electron density, velocity, temperature, and the electrostatic 
potential, respectively. D(u) = {dij)fj^i, dij = ^{diUj + djUi). The constants 
and /i are the viscosity coefficients with p, > and 2/i -I- Niy > 0. cy > is the 
specific heat constant, k > the heat conductivity coefficient, and A > the scaled 
Debye length. The pressure function P{p, 0) takes the form 

P{p,9)^Rp9, R>0. (1.5) 

Without loss of generality, we assume cy = i? = 1 for notational simplicity. The 
Navier-Stokes-Poisson system (|l.ip - (|1.4l) can be used to describe the dynamics of 
plasma, where the compressible fluid of electron interacts with its own electric field 
against a charged ion background, see Degond [3]. 

The purpose of the present paper is to investigate the quasineutral limit of the 
compressible Navier-Stokes-Poisson system (|l.l[) - (|1.4p . We shall prove rigorously 
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that, as the Debye length A ^ 0, the solution of the compressible Navier-Stokes- 
Poisson system converges strongly to the strong solution of the incompressible 
Navier-Stokes equations plus a term of fast singular oscillating gradient vector fields 
as long as the strong solution of the latter exists. Moreover, we also consider the 
convergence of the compressible Navier-Stokes-Poisson system (|l.ip - (|1.4l) to the in- 
compressible Euler equations by performing the combined quasineutral, vanishing 
viscosity and vanishing heat conductivity limit, i.e. A — > and fi, v, k —>■ 0. 

We first give some formal analysis. We use the subscript A to indicate that the 
unknowns are dependent on A and set 4>x — X^\. Thus, we can rewrite the system 
(inil-dlll) as 

dtpx + div(pAUA) = 0, (1.6) 

px{dtux + (ua • V)ua} + V(pa6'a) + jPx^cpx - A^Aua + {i^ + /i)VdivuA, (1.7) 

PxidtOx + (ua • V)6Ia} + PA^'AdivuA 

= kA9x + t/(divuA)2 + 2/iD(uA) : 0(ua), (1.8) 

- AA0A =Px-l- (1.9) 

The system (|1.6p - (jl.9|) is equipped with the initial data 

px{x,0) = pox{x), UA(a;,0) = UoA(a;), 6'A(a;, 0) = 6'oA(a;)- (1-10) 

Letting A — > formally in the Poisson equation p.9p . we have pA = 1- Moreover, 
if we assume that 

UA ^ V, 6x^9 

as A ^ 0, we may expect that the compressible Navier-Stokcs-Poisson system ()1.6|) - 
(|1.9|) converges to the incompressible Navier-Stokes equations (see [17]) 

V-v^O, 

dtv + (v • V)v + Vn = /iAv, 

N (1-11) 

dt9 + (v ■ V)9 ^kA9+!^Y1 + ^J^')"' 

h3 = l 

as the Debye length goes to zero, where VII is expected to be taken as the limit of 
the singular electric field and the gradient of pressure together. Furthermore, if we 
let /Lt — > and /c — > in (jl.lip . it yields the incompressible Euler equations 

V-v==0, 

5tv + (v V)v + vn = 0, (1.12) 
dt9 + (v • \7)9 0. 

Recently, there are many progresses on the quasineutral limit of the compress- 
ible isentropic Navier-Stokes-Poisson system (i.e. the system (ll.6p . (jl.7p and (|1.9p 
with the pressure Pa — o.p1,-f > l,a > 0), Wang [23] studied the quasineutral 
limit for the smooth solution with well-prepared initial data. Wang and Jiang [24] 
studied the combined quasineutral and inviscid limit of the compressible Navier- 
Stokes-Poisson system for weak solution and obtained the convergence of Navier- 
Stokes-Poisson system to the incompressible Euler equations with general initial 
data. In [24], the vanishing of viscosity coefficients was required in order to take 
the quasineutral limit and no convergence rate was derived therein. Ju, Li and 
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Wang [11] improved the arguments in [24] and obtained the convergence rate. Do- 
natelli and Marcati [4] investigated the quasineutral limit of the isentropic Navier- 
Stokes-Poisson system in the whole space R'^ and obtained the convergence of weak 
solution of the Navier-Stokes-Poisson system to the weak solution of the incom- 
pressible Navier-Stokes equations by means of dispersive estimates of Strichartz's 
type under the assumption that the Mach number is related to the Debye length. 
Notice that their arguments can not be applied to the periodic case since the dis- 
persive phenomenon disappears in this situation. Ju, Li and Wang [10] studied the 
quasineutral limit of the isentropic Navier-Stokes-Poisson system both in the whole 
space and in the torus without restriction on the viscosity coefficients. 

However, there is no analysis on the quasineutral limit of the compressible non- 
isentropic Navier-Stokes-Poisson system yet. In the present paper, we shall consider 
the general ill-prepared initial data for the system ()1.6p - ()1.9|) . so the fast oscillating 
singular term will be produced by the non-divergence free part of initial momentum, 
and has to be described carefully in order to pass into the quasineutral limit. 

In order to describe the oscillations in time, we introduce the following group 
C — e'^^,r G R, where L is the operator defined on the space H — (L^(T^))^ x 
{V^,V^ e H^T^)} by 







L[ 1=0, if divw = 0, 



(1.13) 



Then it is easy to check that e^^ is an isometry on space H^{T^) x H'^{T^). Let 
us consider the evolution of velocity and electric field. From (|1.7[) and (|1.9p . it is 
easy to obtain the following equation 

9tV0A-^QuA = -Q(uAV-(V(/.A)), (1.14) 

where the operator Qv = VA^^V • v is the Leray's projector on the space of 
gradient of vector field v e (L^(T^))^, which is defined as follows 

Qv = VA-^V • V, T'v = (/ - Q)v, V ■Pv = 0. 

We project the momentum equation (|1.7[) on the "gradient vector fields" to obtain 

dtQu^ + ^V0A = - Q((UA • V)ua) 

- Q( — VPa) + mQ(Aua) + + M)Q(VdivuA) 

PA 

+ /iQf(— - 1)Aua) +(!. + /i)Q((— -l)VdivuA). (1.15) 



Define 



Then we can rewrite the system (|1.14l) - (|1.15l) as 

x) [ fci 



d^V>^ = ^{-T)( J" (1-16) 
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with 



ko = - Q((uA • V)ua) - Q(— VPa) + mQ(Aua) + (i^ + M)Q(VdivuA) 

PA 



+ mS((^ -1)Aua) +(i^ + m)q((^ -l)VdivuA), (1.17) 

fci =Q(uaV-(V0a)). (1.18) 

Now we can construct the oscillating terms as follows. Let v e C([0, T]; H'^{T^)) 
be a divergence free function. Consider the following linear system 

dtVq + (v ■ V)Vq + (Vg • V)v + vAq) - {^l + J//2) Vdiv(Vg) = 0, 

^ ^ ^ (1.19) 

dtVp + 2 • ^)^P + (^P • ^)"^ + ^^p) - (P + i^/2)Vdiv(Vp) = 

with initial data 

{\/q{x,0),\/p{x,0)) = (Quo(x),V0o(a;)). 
It is direct to prove that there exists a unique global smooth solution (Vg, Vp) to 
the oscillating system ()1.19p satisfying 

||(Vg, Vp)(i)||H.(T") < C(r)||(Quo, V0o)||^^.(T"), (1.20) 

where C(T) > is a constant depending only on T. 
Define 

Uosc(a;,<) \ _ / Vg(x,t) \ , , 

Before stating our results rigorously, we first recall the local well-posedness result 
on the initial value problem for the incompressible Navier-Stokes system ()l.lip in 
multi-dimension. One can refer to [17] for the proof. 

Proposition 1.1. Assume that s > N/2 + 1 and 

v(a;,0) = vo(a;) e divvo = 0, 

9(x,0) = Ooix) e inf 9o(x) > 0. '^^'^^^ 

Then there exists some time T*{0 < T* < +oo) such that the initial problem (|l.lip 
and (jl.22p admits a unique strong solution (v, 0) satisfying, for any T < T* , 

veC^([0,T],i7^+3-), z-0,1, ||v(t)||H.+3 <Co||vo||ffe+3, (1.23) 

eeC\[0,T],H'+^~'), i = 0,l, \\9{t)\\Hs+s<Co\Wo\\Hs+^ (1.24) 

with Co > a constant. Moreover, if N = 2, the initial problem (jl.lip and (ll.22[) 
admits a global unique strong solution (v, 9) € C*([0, oo), H^^^^^), i = 0,1. 

Our main results of this paper read as follows. 

Theorem 1.2. Let < T < T* defined in Provosition \ 1 . 1\ and suppose that (v, 9) € 
C"([0, T], H'^^'^^^), i = 0, 1, s > N/2 + 2, be the unique strong solution of the initial 
problem (|l.lip and (jl.22p . Assume that the initial data (pqx(x),Uqx{x),9q\{x)) 
satisfies 

Poa(x) = 1-AA0oa(x), inf poA(a;) > 0, V(box e H'+^T""), (1.25) 
UoA e H%T^), 9ox{x) e i/^(T^), inf 9oxix) > 0, (1.26) 
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and 

llT'uoA - vollff. + IIQuoA - Quollff. < CX, (1.27) 

Wpoxix) - 1 + XAMx)\\h^ < CX'', pox - OoWh^ < CX (1.28) 

for some constant C > 0, where 0o and Uq are defined by (j2.ip . Then there is 
a small constant 6t > such that, for any X G {0,5t], the initial value prob- 
lem for Navier-Stokes-Poisson system (jl.6p - (|1.9p admits a unique classical solution 
(px,ux,9x,'p\) on [0,T] satisfying 

sup \\{px,ux,ex){t)\\Hs + sup II V0A(t) II (1.29) 

0<t<T 0<t<T 

uniformly with respect to X. Moreover, it holds that 

sup {wipx - imu^ + ii(uA - V - uosc)(oiih= + widx - omwn^} 

0<t<T 

+ sup ||(V0A- V0osc)(i)||He+i < C2A (1.30) 

0<t<T 

with C2 > independent of X. 

If we further perform the combined quasineutral, vanishing viscosity and vanish- 
ing heat conductivity hmit, i.e. A ^ and p,i>, k 0, we obtain the convergence of 
the Navier-Stokes-Poisson system ()1.1|) - (|1.4|) to the incompressible Euler equations 
(frT2| . Namely, 

Theorem 1.3. Let < T < T** and suppose that (v, 6*) G C'{[0,T], H'+^-'^), 
i ^ 0,1, s > N/2 + 2, be the unique strong solution of the initial problem (|1.12p and 
()1.22p . where T** is the maximal existing time of{-v, 9). Assume that the initial data 
{pf)x{x),'VLQx{x),Oox{x)) satisfies the conditions (|1.25p - ()1.28p . Then, there is a small 
constant 5t > such that, for any X £ (0, 5t\, the initial value problem for Navier- 
Stokes-Poisson system p.6p - (|1.9p admits a unique classical solution (px, Ux, Ox, 4>x) 
on [0, T] satisfying 

sup \\{px,nx,9x){t)\\H^ + sup \\y(^x{t)\\H^+i < C3 (1.31) 

0<t<T 0<t<T 

uniformly with respect to X as fj,,^, k 0. Moreover, it holds that 

sup {wipx - imu^ + ii(uA - V - uosc)(oik= + m - ornwH^} 

0<t<T 

+ sup ||(V0A- V0osc)(i)||He+i <C4A (1.32) 
o<t<r 

with C4 > independent of X. Here (v, 9) is the unique strong solution of the initial 
problem (|1.12p and (|1.22p . and (uosc^osc) is the fast singular oscillating gradient 
velocity vector field and electric field defined by (I1.19P and (ll.2ip with p — v = 0. 

Remark 1.1. The method developed in this paper can be applied to the situation 
when the doping function is a perturbation of a constant state 

C{x) = l + Xg{x) 

with g{x) e C^(T^), a given function, satisfying J^jy gdx = 0. 

Remark 1.2. We believe that the method developed in this paper can be also applied 
to investigate the quasineutral limit problem to more complex model such as the 
full Navier-Stokes-Poisson system with more general pressure, which will be studied 
in a forthcoming paper. 
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The proofs of Theorems 11.21 and 11.31 mainly consist of three steps. First, we 
apply the homogenization technique to construct the approximate solution to the 
classical solution (if exists) of the system (|1.6p - (|1.9p . Then by using the theories of 
symmetric quasilinear hyperbolic system and the estimates of second order elliptic 
equations, we show that the remainder term exists in the same time interval as 
the approximate term for fixed small A > 0. Moreover, we obtain the uniform 
estimates with respect to A (the uniform estimates with respect to /i, v and k can 
also be obtained by further analysis). These facts are sufficient for us to complete 
the proofs of Theorems 11.21 and 11.31 

It should be noted that the quaineutral limit is a well-known challenging and 
modelling problem in fluid dynamics and kinetic models for semiconductors and 
plasmas. In both cases there exist only partial results. In particular, the quasineu- 
tral limit has been performed in Vlasov-Poisson system by Brenier [1] , Grenicr [5] , 
and Masmoudi [18], in Vlasov-Poisson- Fokker-Planck system by Hsiao, Li and 
Wang [7,8], in Schrodinger-Poisson system by Puel [21], Jiingel and Wang [13], 
and Ju et al. [9], in drift-diffusion-Poisson system by Gasser et al. [6], Jiingel 
and Peng [12], Wang et al. [25]. For the hydrodynamic model, besides the re- 
sults mentioned above for the Navier-Stokes-Poisson system, there are also many 
results on Euler-Poisson system, for example, for the isentropic Euler-Poisson sys- 
tem [2,19,22,23] and for non-isentropic Euler-Poisson system [16,20]. Li and Lin [14] 
considered the quasincutral limit to the isentropic quantum hydrodynamical model 
with the help of modulated energy method for general initial data. 

Before ending this section, we recall the following Moser-type calculus inequali- 
ties which will be used frequently in the sequel. 

Proposition 1.4 ([15] Moser-type inequalities). (1) For f,g G n and 

\a\ < s, it holds that 

II^"(/.9)I|l^ < CML^WD'gh^ + Ml^WD^Wl^). (1.33) 

(2) For f e H'',Df e L°°,g e i?""^ n L°° and \a\ < s, it holds that 

\\D"ifg) - fD"ig)U2 < CsiWDfh^WD^-'gU^ + Ml^WD^ /Wl^). (1.34) 

Notations. In this paper, C and C'i{i = 1,2,...) denote the generic positive 
constants, which may change from line to line and are independent of A. C(T) 
and Ci{T) denote the constant depending on the time T. denotes the standard 
Sobolev space iy*'^(T^). For the multi-index a — (ai, . . . , um), we denote 13" = 
d^l---d^^: and|a| = |ai| + ----f |a^|. 

The rest of this paper is arranged as follows. In Section 2, we construct the ap- 
proximate solutions to the problem (|1.6p - (|1.10p . In Section 3, we establish the local 
existence of solution to the remainder system and obtain the uniform estimates. 
The proofs of our main results are given in Section 4. 



2. Construction of approximate solutions 

In this section we shall construct the approximation to the system (jl.6|) - (|1.9|) . 
Noticing the fast singular oscillating vector fields (uqsc, V(/)osc) obtained by (|1.2ip . 
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we find that the fast smgular osciUating vector fields (uosc, V(/)osc) satisfy 
( 1 / 

(9tUosc + ^Q{{v ■ V)Uosc + (Uosc • V)v + VV • Uo 



(/^ + i//2)Vdivu, 



osc r -^V0osc — 0, 



5tV</)osc + 2 2( (V ■ V)V0osc + (V0OSC • V)v + vA0o 



(2.1) 



(/i + t//2)VA0o 



0, 



(UqSC (-^; 

Thus it is natural to define 

Pose — A(^osc- 



We conclude that the fast oscillating part (posc, Uqsc, 0osc 
initial value problem 

1 



) satisfies the following 



dtPosc + [V + Uosc] • Vposc + + ^Posc)V • Uosc = fc2, 

9tUosc + ([v + Uosc] • V)Uosc + (Uosc ' V)v + jVcposc = ^3, (2.2) 
^0osc Pose; 

. Posc(a:, 0) = -A(/)o(x), Uosc(2:, 0) = Quo{x), 



where 



V • (posc[v + Uosc]) + 2 V • ((V • V)V0osc + (V0OSC ' V)v + vA^osc) 
- (M + i^/2)A20osc, 

-Q((V • V)Uosc + (Uosc • V)V - vV 

'^osc) (U-osc • V)Uosc 



(2.3) 
(2.4) 



+ V{[W ■ V)Uosc + (Uosc • V)v) + + I//2)VdivUosc- 

Moreover, by virtue of (|1.20p and (|1.2ip . we obtain that 

l|fc2||i/=-2(T") + ||fc3||i/=-2(T") < C\\{Vcj)o, Quo, Vo)||h=(T")> (2-5) 

where the constant C > is independent of A. To approximate the classical solution 
W = {p\,u\,9\,(j>\)'^ of the initial value problem (|1.6p - (|1.10p for small A, we still 
need to introduce an additional correction term 

^^cor ( Apcor 1 Ucor 7 ^cor 1 f/^cor ) 

By utilizing the fast singular oscillating part and the given functions /c2 and fcs, 
we can construct (pcor, Ucor, 6'cor, </'cor) by solving the following linear initial value 
problem 

' l9rUcor + V0cor = h, 



^rVt/icor - Ucor 

Pcor A0CO1'; 

drf^cor — ^5 J 
_ (Ucor,V(^cor,0cor)(a;,O) = (0,0,0), 



(2.6) 
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where 



ki = - -V6 + i-J. Auosc + {p + t^)VdivUosc, 

fc5 = - Uosc -S/e -0V ■ Uosc + I^(divUosc)^ 



N 



Here we recall that (v, 9) is the solution to the system (II. lip . 

By virtue of (fL^ . dOSD, PTM)) and it is easy to prove the following 

existence results of solutions to the problems (|2.2p and (|2.6p . 

Proposition 2.1. Let T > Q,T < T* be given. Letv,e £ C"([0, T], i = 

0,1, s > 1 + N/2, be the solution to the initial value problem p. lip and (|1.22p . 
T/ien i/ie problem (|2.2p admits a unique classical solution (posc, Uqsc, V^osc)"'" /o?^ 
i e [0; r] satisfying 

||Posc(i)|l_ff-' + 2 + ||(Uosc, V(/)osc)(0lli/=+3(T«) < C't, (2.7) 

and i/ie problem (|2.6p admits a unique classical solution (pcor, Ucor, fi'cor, V(/)cor)"'" 
/or i G [0, T] satisfying 

V0cor)(T)||H=i+2(T«) < Ct, (2.8) 

where Ct > depends only on T and the initial data (vq, 6'o, Quq, V0o); but is 
independent of X. 

According to Propositions 11.11 and 12. 1[ we can make the following asymptotic 
expansions of the solution (px, Ux,0Xt4'\) 

px{x, t) = \ + Xposc{x,t) + A^(An(a;, t) + pcoi{x, t/X)) + X'^picmix, t), 

Ux{x,t) = V + Uosc{x,t) + XUcor{x,t/X) + XUrcm{x,t), 

6x(x, t) — d{x, t) + X9coi{x, t/X) + X9r 

cm ^) •) 

(j)x{x,t) = cj)osc{x,t) + A(n(a;,t) + cl)cor{x,t/X)) + X(f>rcm{x,t). 

Substituting ([^ into the Navier-Stokes-Poisson system prB)) - P3)) . using (fLTTjl . 
(j2.2p and (j2.6p . and by tedious but direct computations, we can show that (prcmj 
Uromj ^romj '^rom) solves the foUowing initial value problem 

dtPrcm + Ua • Vpiom + ^PAdivUiom = ^0, 



rem ~ (uA • V)U 

rem 

+ A— Vp 

rem rem 

p\ 

- /iAUiom - (M + J')VdivUrcm = ~l-'^'f'i-cm + fo, ^, 

(2-10) 

C^t6'rem + Ua • V^rcm + 6'AdivUi.cm ~ KA9^cm = Al'(divUrem)^ 
fiX ^ 



2 



- A0rom = Pr. 
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with initial data 

Prem(a;, = ^ [pox{x) " 1 + XA(j)o{x)] - AU{x, 0), 
Urom(a;, 0) = ^ [uoA(a;) - vo(x) - Quo{x)] , (2.11) 

^rom(a;,0) \ [eox{x)^0a{x)]. 
A 

In (|2?T0| , we denote 

ho = - Urom • Vposc - PromV • (Uqsc + AUcor) - V ■ (poscUcor) 

- (v + Uosc + AUcor + Au 

- AHt - (V(An))(v + Uosc + AUcor + AUrcm) 

- Andiv(Uosc + AUcor), (2.12) 

fo=foi+fo2, (2.13) 
go =901 +.902 (2.14) 

with 



foi = - ((Ucor + Urcm) ' V) (v + U^sc) - ((v + U^sc + AUcor + AUrem) ' V)Ucor 

- — V(posc + A(An + peor)) ^ V^cor, 
PA 

f02 =pAUcor + (p + i^)VdivUcor 

\Posc 

+ A(An + pcor) + Aprcm) A(V + 

Uosc + AUcor + AUicm) 

PA 

(Pose 

+ A(An + Peer) + Aprom) Vdiv(v + 

Uosc + AUcor + AUrcm): 

PA 

gOl = - (Ucor + U,cm)V6' - (v + Uqsc + AUcor + Au,cm)V6'cor 

- (^'coI■ + 6lrcm)divUosc + ^AdivUcor, 

g02 =KA6'cor - — (Posc + A(An + Poor) + Api.cm)^(^ + A^cor + A6'rcm) 
PA 

+ 2l^divUosc (divUcor + divUrem) + Al^(divUcor)^ + 2Ai^divUcoi divUi.cm 

N 

+ p ^ {^^VJ + d.jv^ + diul^^ + djul^^){d^ul^, + djul^, + di 
X ^ . ^ . 



vPosc 

+ AAn + Apcor + Ap 

PA 



N 

Kdivu,)2 + | ^(a,ul + 9,ul)^ 



If we denote 

f^rcm • (Prcmi U] 



r, \T 

rem ; ^rcm ) : 
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the problem p.l0p - ()2.1ip can be rewritten as follows 

N 

dtUrcm + ^ Aj {x, t, Urcm)dxj Uicm - M^Uiom - (M + i^) VdivUrcm 

J = l 

- KA^rem = Al^J + +\b + F{x, i, [/rem), (2.15) 

^4^rcm Prcirn 

Urcm{x,0) = (prom(a;,0),U,.cm(2:,0),6',.em(2:, 0))"^ := Urcma{x). 

Here the matrices Aj{j — 1 , . . . , N) is defined as 



I yxej, 

MieJ O 
\ exe, 



and 

Urcm — (O5 Urcm, 0) 7 ^rcm — (O5 ■ ■ ■ 7 ^rcm) 5 

J= (0,...,0,(divu„,„)')^, F= (/io,fo,5o)^, 

. N . T 

G = I 0, . . . , 0, ^ (a,<,„ + 9,<em)' J , S = (0, -V(^rem, 0)^. 



3. Local existence of solution to the remainder system (|2.15p 

In this section we study the local existence of smooth solution to the remainder 
system (|2.15p . our result reads 

Theorem 3.1. Let T > 0,T < T* be given and v,0 e C'{[0,T], H'+^~'),i = 
0,1, s > 2 + N/2, be the solution to the problem (jl.lip and (|1.22p . Then there 
exists a constant St > such that for any X £ {0,St], the initial value problem 
(|2.15p admits a unique classical solution (C/rom, '/'rem) in [0,T] satisfying 

sup {\\{Xpi- 

cm, Urcm, ^rcm) (0 II -f^^ ^~ II ^'/'rcm (0 1 1 + i ) — ^(-^), (^■-^) 

0<t<T 

where C{T) is a positive constant independent of X. 

The proof of Theorem [STT] proceeds via a priori energy estimates and the classical 
iteration scheme. The crucial step is to show the following energy estimates which 
can be obtained by performing the refined energy estimates for the quasilinear 
symmetric hyperbolic-parabolic system and the Poisson equation. 

Lemma 3.2. Let T > be given and s > N/2 + 2. There exist positive constants 
5t,M,M such that the classical solutions (f/rem, 0rem) to the initial value problem 
(|2.15p satisfies 

sup (||(A/9i.om,Ui.cm,^'rom)(i)||lfa + || V(/)rcm (0 II + i ) 
0<t<T 
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sup 

0<t<T 



(||A9tPrem(i)||ff'-i + || A9t U,.cm (i) 1 1 + \\dtO,cm{t)\\H''-^ 

+ \\XdtyAemmHA <M (3.3) 



uniformly with respect to A G (0, Sj 



Proof of Lemma \3.SX We assume a priori that the classical solution to initial value 
problem ()2.15p satisfies (|3.2p and (j3.3p . Then our task is to determine these un- 
known constants by energy estimates. 

Noticing the matrices Aj(x, t, t/rcm), j ~ \, . . . ,N can be symmetrized by 



O 

O pxInkn O 
o ^ 



Ao{x,t, C/rem) = 

we rewrite the system (|2.15|) in the following form 



N 



cm rc 

i=i 

- + I^)pAVdivUi.om - -^A^rem 

PA 

= Ai. J + +\b + F{x, t, c/rem), 

Z A 

^^rcm — Prcmi 
^ [/rem (2^,0) = C/i.emo(a;), 

where ylj = AoAj,j = 1, . . . , N are symmetric matrices given by 

/ XexCj 
AjiUrcni) ^u{Aa{U,cnd + [ ^^^xej O pxej 











and 



J:- AoJ- (0,...,0,^(divUrem)' 
^ fA 

. N s 

G := AoG = f 0, . . . , 0, g- ^ (5.<m + djul^nf) 

1,1=1 



Pa t/A 



F /IF f ^^xhp Pxgo\ 

F := = I ,Pxto, 1 

\ Px »x ' 



(3.4) 



Next we perform energy estimates for the classical solution to the system (|2.15p 
with initial data (|2.1ip . Define the canonical energy by 



If/. 



rem 1 1 • ~ 



(Ao(t/rcm)^rcmi Urcm)dx. 
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Multiplying ()3.4p i by J7rcm and integrating the result by parts, we get the basic 
energy equality of Friedrich's 



^||t/'rem||l; 



2k I ^\W9rcn,\^dx 



y(rC/rom, U,-cm)dx + 2\v j ^ (divUrom)^6'rcmC?a; 

N 

+ 2 [(AoF,Urom)dx + Ri, (3.5) 



where 



Rl =2{fi + v) j {px - l)VdivUromUrcmda; + 2^ y (pa - l)AUi.cmUromrfa; 

-2k /v(^)V0,.em0,cmdx (3.6) 



and 

r=.(at,v)-(Ao,A,...,^3). 

Since /i > 0, 2/i + A^i^ > 0, there exists a positive constant such that 

H j \V\^.,^^\^dx + {^l + v) j Idivurcml^dx > Ci j \Vurcm?dx (3.7) 

in view of /(divUicin)^rfa; < / \^ '^rcm\'^ dx . Notice the fact that there is a > 
such that for A £ (0, A^] it holds that 

< P- < 1 + Aposc + A^An + AVcor + AVrcm < P+, (3.8) 

Q<e^<0 + xe,or + A^^,.cm < o+, (3.9) 

where p± and 9± are positive constants. Thus, the matrices Aq and Aj,j — 
1,...,N, together with their derivatives are continuous and bounded uniformly. 
Moreover, Aq is uniformly positive definite, i.e. there exists a constant co > such 
that 

{Ao{U,^m)Uro,n, f/rcm) > Co(AV?cm + U^cm + Cm) (3-10) 

for all U^cm- 

Now we estimate the terms on the right-hand side of (|3.5p . Since F is bounded 
there exists a generic constant A/q, independent of (prom, Urcm, ^'rcm, 0rcm) and A > 
0, such that 



(Ff/,„„, C/rcm)dx < Mo(l + A(M + M))\\U,.,n.\\E. (3.11) 
By Sobolev's embedding inequality and the inequality (|3.9p we obtain that 
f 1 ^ f I - 

2Xv j —{divU,.cm f9rcmdx + Xp ^ J — ((9iM^p„i + 9-,uJ„-,i)^6'remda: 
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<\MoM{2fi +U) j (IVUremP + l^rcmHrfa;. (3.12) 

By integrating by parts, Cauchy's inequality and the equation for prcm in (|2.15p . 
the forth term on the right-hand side of (|3.5p is estimated as follows 



rem '-^rcm 



= - 2 y dtPrcm4'i-cmdx - 2 J + "osc + AUcor + AUrcm) Vprcm<^rcmda: 

+ 2 J J /lo(/)romrfa; + 2 y V(posc + A(An + Pcor) + Xprcm)Urcm4>rcmdx 
< - 5t||V(/.rem||i2 + A/o(l + AAf ) (|| V0rcm 1 1 ^2 + ||i7rcm|||) + £1 1| Vu^om || ^2 (3.13) 

for some sufficiently small constant ei > 0. 

Now we deal with the term By integrating by parts and using Sobolev's 
inequality, we get 

2{p + '^) J (pa - l)VdivUremUromrfa; ~^ '^^ j ^ I) ^U^ciaVi^cmdx 

<XMo{M + l){2fi + i^) I dVurcml' + |u„,„|')rfa;. (3.14) 



In view of p.8p . (|3.9p and Cauchy's inequality, we obtain that 

-2k y" V^V6'rcm^'remda; 
= - 2k J ^^^Ve^-cnAcmdx + 2k J ^^^^^ V6xS7e^.cnAcmdx 



<XMo{M + 1)k J (|V^?rcmP + \Orem\^)dx 

+ MoK J \9remfdx + e2K J \Vercmfdx (3.15) 

for some sufficiently small constant £2 > 0. 

The estimate of the fifth term on the right-hand side of (|3.5p is tedious but 
straightforward. In view of the definitions of Hq, fo, and go in (|2.12p - (|2.14p . and the 
Propositions 11.11 and 12.11 we get 

2A^ [ —hoPrcmdx + 2 [ pxtoiUrcmdx + 2 [ ^goi9rcmdx 



P\ J J ox 

and 



<Mo\\U,,^\\l + Mo (3.16) 



2 / pxfo2 ■ Uromrfa; + 2 / ^go20rcmdx 



ex- 

<X{2fi + J/ + k)Mo(1 + M) J (|u,em|^ + IVUrcml^ + \Ve,^n,f)dx 

+ {2p + u + k+l)AIo. (3.17) 
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We choose St sufficiently small such that, for A G {0,5t], 



\Mo{M + l)(2/i + I. + «) < min {|-, -g-} := 771. (3.18) 

Choosing ci and €2 sufficiently small and combining p.7p - p.l8p with p.Sp . we 
obtain that 



<Mo(l + \{M + A?))(||[/,.en,||| + ||V(^rem|li2) + St? J i\u,,^\^ + \e,,^,\^)dx 

+ kMo / |6lremprfa; + (2/i + i/ + K + l)Mo. (3.19) 



Next we shall obtain the energy estimates of higher order derivatives for the 
classical solutions to the initial value problem (|2.15p . For the multi- index a with 
1 < \a\ < s, we take the operator _D" to (|2.15[) and multiply the resulting equations 
by Ao to obtain 

N 

i=i 



^{p + iy)pAVdivi?"Urcm - ^AD"0~,em 

PA 

= Ai.Ao(C/,em)^"^ + ^Ao(C/rcm)i?"G + ^ (C/,cm)i?"B 



(3.20) 



+ Ao(C/„m)^"i^ + i?", 
- A£)"0rcm = D°'Picm 

with initial data 

7^"C/rcm(a;, 0) = I?"{/rcmo(a;), (3.21) 
where i?" consists of the commutating terms as 

N 

H" = -^Ao(t/,em)(^"(A,(t/rem)5,^t/rem) - A,([/,em)a,^^"f/,em). 

i=i 

Taking the inner product between (|3.20p i and D^Uicm, we have the following 
differential equality 



-I 



i)"C/rem(i)|||; + 2M j \V D'' Vi,,^\^ dx + 2[^1 + u) J \div D" U,,^\'' dx 



2n I ^ 



2 



dx 



: y (ri?"C/,c,„, D''U,,,,,)dx + 2Az. J{Ao{U,,^)D"J, D"U,,,n)dx 

+ A/i y" (Ao(C/rom)i?"G,I?"C/,.cm)d.T + | ^ ( Ao (iJrcm)^"^, i?" iJrem)^^; 

+ 2 / (Ao({/rcm)i^"i^(t),i5"C/rem)da; + 2 / {H'^ (t) , D"U,em)dx + i?2, (3.22) 
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where 



+ 2(/i + 1^) J (px - l)VdivD"u,e,„i:'"u,cmda;. 
It is easy to see that we also have the foUowing estimate 

M J \VD°'u,^,^\^dx+{fi + i') J |div£)"urem|^rfa; > |VD"Urempda; (3.23) 

for some constant ^2 > 0. 

Now we deal with the right-hand side of (|3.22p . In the following the generic 
constant Mq may depend on T and s. By integrating by part, Sobolev's inequality 
and Cauchy's inequality it holds, similar to (|3.1ip and (|3.14p - p.l5p . that 



{TD^U,,^, D"U,,m)dx < Afo(l + A(M + M)) p"C/,e„,|||5 (3.24) 



and 



-^rcm I 



R2 <XMa{M + l)i2^i + ly + J (|VD"u,„„P + |i5"+i0,.cmP + |i?"u, 

+ |i?"0rcm|')rfx + AfoK J \D'' 6,,^]'' dx + 6kJ |i?"+l0rom|' (3.25) 

for some sufficiently small constant S > 0. 

By the definition of Ao,G and J, it follows from the Sobolev's inequality that 

Xiy J{Ao{U,,,^)D''J,D"Urcn.)dx + 2Xfi J {Ao{Urcn.)D''G, D''U,,^)dx 

<XMo{2fi + V) f ||(divU,e„02||^„ + II V (C)^<„, + \ ll^rcmlU" 

< A Afo (2/1 + I/) 1 1 U,em 1 1 1= 1 1 Urcm | ] |„ + 1 1 1 B,^,^ \ \ ffc 

<AA/o(2/i + t^)||u,cm||4„||u,em|||„ + l||ercm|||c + AMq (2/t + /^) |1 U,.em || + i 
<AA/A//o(2/t+ J^)||u,e,„|||„ + i||0rcm||ff= + AA/o(2/i + J/)||Urem||ff. + l. (3.26) 

We deal with the fourth term on the right-hand side of p.22p . From (|2.15p . we 
can easily get the equation for ZJ^prcm, 

fti?"prcm + UA • Vi?"prcm + ^PAdiv/^^U^cm = i?"/lO + /l" (3.27) 



with 



D«(ua • Vpren,) + Ua • VD"p,e„, - ^Z?" (pAdivU,.„„) 
ipAdivD^Ui-om- 



In view of (|3.27p and the Poisson equation (|3.20p 9. we get 

^ y" (Ao(C/rcm)i?"S,i?"C/rcm)da; = -f / PAVi?"0rc,„i^"u„,„dx 
= ^ / PAdivi:i"u,.om-D"0rcmrfa; + ^ _/ VpA-D^Urom-D^^rcmC^a; 
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+ 2 J D"hoD''(j),^n,dx + 2 J V(posc+A(An + Pcor) + Aprcm)£'"u„m£'"</'rcmrfa: 
+ 2 h'^D^chrrvndx 



0<\P\<\a\ 



2 

cm 1 1 E 



+ 63 J IIVD^Urcmlrda; (3.28) 

for some sufficiently small constant £3 > 0. 

The fifth term on the right-hand side of p.22p is very tedious. The main tech- 
niques involved are Leibniz's formula, Moser-type calculus inequalities (|1.33p - (|1.34p . 
and Sobolev's embedding inequalities. Actually, after the tedious computations, we 
finally obtain the following estimate 



2 J {AoiU,,m)D''F{x,t,Urcra),D''Urera)dx 

<A(2/i + + k)Mo(1 + M) 



J2 {\\VD''u,,^,\\l, + \\\7D^9,,^,\\l,) 

0<\P\<\a\ 



,|2 
rem 1 1 E 



0<\l3\<\a\ 

The commutating term H" can be bounded by 

;i/"(t),i?";7,em)dx 



+ (2/i + zy + K + 1)Mo. (3.29) 



< E Mo(l + AA/)||I?'^[/re„,|||-f ||I?"t/,em||| + Afo. (3.30) 

l<\0\<\a\ 

We now re-choose 6t sufficiently small such that, for A G (0, St], 

AsAfo(A/ + l)(2A* + !^ + «;) <min{^,^} :=r;2. (3.31) 

Let 

Ht) = X'WprcraWh + l|Urcm||ff. + ll^rcmUlf^. (3.32) 

Taking S and £3 small enough and combining the estimates (|3.24p - (|3.30p with p.22p 
and (|3.19p . we obtain that 

{a/o (a/o(1 + A(Af -f Af )) -f 37/ + A/oK + A(2/i + t/)A'/A//o||u,.en,|||,+i) 

X (co$(r) + l|V0,emlllf.(r)) ]dr + co$(0) + || V(/.,em(0)|||j. + A/o(2Ai + u + k)T, 

(3.33) 
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where ^ = niinj^i, ^2} and 77 — max{?7i, 772}- By virtue of Gronwall's inequality, we 
obtain that 



Co$(i) + ||V0rem||H= < (co$(0) + || V^^em (0) 1 1 l^a + Mo{2^1 + ly + k)T) 

X exp|Mo I A/o(l + A(M + M)) + 877 + 7\foK + A(2^ + i^)MAfo||urom|||,+i rf?-}- 



(3.34) 

From p.2p and Holder's inequality, we have 



A(2/i + v)MMq / ||urcm|l|j.+idr < \Mo{2^i + iy)M^T^. (3.35) 
In view of (|1.27p and p.28p . we obtain that 

A'||Prem(0)l|l^. < CA^ || U,em(0) || 1^. + ||0rem(O) < C (3.36) 

and 

\\^Acm\\h<C. (3.37) 

We choose 6t sufficiently small such that, for A G (0, St], it holds that 

X{M + M) + \{2^i + v)Mi < 1. (3.38) 

Set 

Li = Mo{2Mq + 3ri + Mqk + MqT^'^). 
Substituting ([SSSll-dSSHl) into ([334ll . we obtain that 

Co$(t) + ||V(/),em||l^. <(CO$(0) + ||V0rem(O)||lf. +Mo(2Ai + J/ + K)T)e^i^ 

<(AfoC' + A/o(2Ai + + K)T)e^i^ L3. (3.39) 
In view of (|3.33[) . we get that 

1/ ||u,em|lH.+irf'^ + -^ / ||0||l,.+idr <LiL3r + AfoC' + Afo(2M + 7^ + K)r. 
^ Jo Jo 

(3.40) 

Therefore (|3.2p is proved if we set 



A^2 (i3 + LiigT + AfoC + Afo(2^ + v + k)T) . max | — , 1, |, (3.41) 

L Co 4 Kp- J 

It follows from (|3?20l) that 

sup [XWdtPrcnMWH'-^ + A||atU,.cm(i) ||h^-2 + ||9t6',e„i(i) 
0<t<T ^ 

+ A||9tV0rem(t)||H=) < A4- (3.42) 

with 

M := (Afo(l + 2Af))i/^ (3.43) 
The proof of Lemma [321 is completed. 

□ 
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Proof of Theorem \3.1[ With the a priori estimates (|3.2p and (13. 3|) , we now start 
the proof of Theorem 13. II We first construct the approximate solutions. Define 



^TT7i+l in-\-l\ ( ■ 



n+1 ,,"+1 ^"+l^T 
rem 7 rem 7 rem ' 



inductively as the solution of linear equations 

N 



(/i + Vdivu; 



i+i '^Pa 



ri+1 



(3.44) 



Z A 



with initial data 



C^rem(a;,0) = C/remO (a^) , 



(3.45) 



where 



p^(a;,i) = 1 + \posc{x,t) + A2(An(a;,0 +Peor(.T,</A)) + A^^emla^, *)> 

u^(a;,<) = V + Uosc{x,t) + Aucor(a;, V^) + ^urom(a;,0> 

ei{x, t) = 0(x, t) + A0cor(x, VA) + ACm(a;, t), 

(j)l{x,t) = </>osc(x,t) + A(n(x,t) +(?!)cor(a;,t/A)) + ACm(a;,^), 

<+i = (0, <+„\ 0)T, - A^B{x, t, C/,em) = (0, ~pl\/c^l+lO), 

An 



J" ^oi^(a:, t, C/;j,^) = (o, . . . , 0, ^(div<„j2^ , 
G" := AoG(x, t, C/,"„J = (O, . . . , 0, ^ ^ ((9.<m)" + {dXeJ' 



It is standard to know that the approximate problem p.44p admits a unique 
solution such that 



<<+ieL2(o,r;i/^+i), 
9tpr,+ieG([0,r];ff^-i), 

and satisfies the uniform estimates 



VC+'et^([0,T];iJ^+i), 
0^0+' eL^(0,T;i/'^+i), 
9*<,+iGG([0,T];i/^-2), 



sup (II (Ap: 

0<t<T 



sup 

0<t<T 



n+1 ,.n+l 

rem ' rem ' 



u; 



n+l||2 
rem + 1 



11^?: 



n+l||2 



^.+1 dt < A'P, 



(A2|ia,p:4+i(i)iii^.-r + A2|ia,u:4+i(t)ii?,.-. + imi+^'mi. 



+ X'\\dtWrZHmHs)<M\ 



(3.46) 



(3.47) 
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It is standard to verify that the difference 

(Pi 



^rcm 7 rem i rem ' 

satisfies 



'^rcm ' V^rcm ) \rr 



n+1 
cm 



n n+1 _ n gn+1 

Frcm' rem rem? rem 



9" (^"+^ 

rem' V^rem 



dtp] 



n+1 



n+1 



-,n+l 



/iAu: 
1 



n+1 



Pi 

(/i + i/)Vdivuj 



rem 



;n+i 



n+1 
rem 



Tl-1 



(3.48) 



I / H H \ 



+ fo{x, t, u';+^,p^+J) - fo{x, t, p^^J 
Observing that, for \a\ < s, 



r^)div<e„. 



<Afo 



j"-OI + |i:»"(G"-G"-OI 



E 

|a|-l=|^| + |7|<s-l 



(3.49) 



Then repeating the previous analysis used in the proof of Lemma 13.21 and using 
the interpolation inequalities, we can show that there is a > such that, for any 
A e (0, St] and s' < s, 

sup (||(Apre1n\<ctr\^7etr)(0llH^.' +l|VCln'WII?,.' + i) 
0<t<T 

Jo Jo 

sup {x^dtpi+lml.'-i + x'\\dtu:+l{t)\\l^,_^, + \\dtei+^\t)\\l^,_, 

0<t<T 

+ x'\\dts/^i+imi^,)<c 

\jcze 

^ G L°°{Q,T; H'') nUp{[0,T]; H"'-^) 



for some constant G > 0. Then the Arzela-Ascoli theorem implies that there exists 
a limit vector function 



(Prcm; ^rcm; ^rcm; ^*^rcm 

satisfying (|3.2p - (|3.3p such that 



sup IKp^em' 

o<t<r 



Prcm ; ^r 



n+1 



- U, 



nn+1 
rem , <^rcm 



VC+ -V0,em)(t)||^.'-2 







as n ^ +00 for any A G (0, i5t]. Furthermore, for N/2 — [N/2] < a < 1, we have 
the convergence 



n+lNT 
rem / 



(Prcm; ^rcm; ^rcm? ^ 4^ri 
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in C([0, T];H'' °") by the standard interpolation inequality. Moreover, by Sobolev's 
embedding theorem, we have 

(Prcm,U„m,^rem,0rem)^ ^ C( [0, T] ; i7 ^' ) H ( [0, T] ; iJ ^'-^ ) 

^ C\[0,T] X T^)nC([0,T];C2(T^)) 

for any A £ (0, (St], where we have used the fact s' > N/2 + 2. Then the exis- 
tence of classical solutions to the initial value problem (|2.15p - p.lip is proved. The 
uniqueness of the classical solutions can be proved easily by energy estimates for 
the difference of any two solutions. Thus the proof of Theorem 13. II is finished. 

□ 

4. Proofs of Theorem 11.21 and Theorem 11.31 

Proof of Theorem ] LSI By the asymptotic expansion (|2.9p . Propositions ll.ll and l2.1[ 
the existence and uniqueness of classical solutions to the initial value problem of 
Navier-Stokes-Poisson system (|1.6p - (|1.9l) is proved and the solution satisfies 

sup \\{px,ux,9x){t)\\H- + sup ||V(/>A(i)||i/a+i 

0<t<T 0<t<T 

+ I|ua||l2(0,T;H= + 1) + ||^'a||l2(o^T;H= + 1) < C{T), 

sup (||5t(pA,UA,0A)(t)|k-^ + \\dtVMt)\\H^ + ^) < C(T,A), 
0<t<T 

where C(T) > is a constant independent of A and C{T,X) > is a constant 
dependent on A. Moreover, it is easy to see that, for A £ (0, St], 

sup \\{px - 1,ua - V - Uosc^'a - 9){t)\\H- 

0<t<T 

+ sup ||(V0a-V(^osc)WIIh=+i <c(r)A. 

Q<t<T 

Thus the proof of Theorem 11.21 is finished. □ 

As far as the combined quasineutral, vanishing viscosity and vanishing heat 
conductivity limit is concerned, we can follow the same lines as the proof of Theorem 
11.21 Recalling the uniformly bounded estimates obtained in Lemma 13.21 we are 
able to get the uniform bound with respect to A, /i, v and k for the solutions. Thus 
Theorem 11.31 can be proved similarly with minor modifications of our previous 
arguments. We omit the details here for conciseness. 
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